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Abstract 

We show that axion domain wall is ferromagnetic in the universe with non- 
vanishing baryon number or lepton number. It is caused by protons and 
electrons bounded to the domain wall with their spins aligned. These bound 
states arise due to attractive potentials generated through pseudo- vector cou- 
plings between the fermions and the axion. Using a model of hadronic axions 
we predict existence of a primordial magnetic field with strength 10 -13 Gauss 

at recombination. 
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The axion [|IJ is a Nambu-Goldstone mode associated with Peccei-Quinn symmetry 
which leads to a possible solution of strong CP problem [[| . The axion gains a mass through 
QCD instanton interactions. These interactions give rise to several possible vacuum expec- 
tation values of the axion field. This fact in general leads to axion domain walls |3|,[J in 
cosmology. Namely in the early universe the axion field can take different vacuum expecta- 
tion values in different causally distinct regions when the temperature of the universe cools 
down to below QCD transition temperature. Thus the domain walls appear between these 
regions. However the energy of the domain walls dominates the energy of the universe so 
quickly that the standard scenario of the early universe is invalidated f|||. In order to solve 
the problem it is helpful to find new properties of the domain wall. 

In this paper we reveal a new feature of the axion domain wall: The domain wall is 
ferromagnetic. This property is caused by polarized neucleons and leptons bounded to the 
wall. We will show that the domain wall produces a potential in which only the particles 
polarized perpendicular to the wall are bounded. The potential has the width of m~ l and 
the depth of m a approximately where m a is the axion mass. Thus the binding energy is the 
order of the axion mass ( 10~ 3 eV~ eV ). As we consider a period when the temperature 
of the universe is higher than that of neucleosynthesis, but lower than the temperature ( 
150 MeV~ 100 MeV ) of QCD transition, relevant particles to the property of the wall 
are protons, electrons, and their anti-particles. These particles form bound states in the 
potential where the spins of the states are aligned in the identical direction perpendicular 
to the wall. Thus as far as the baryon number ( the lepton number ) of the universe is 
non-zero, the wall possesses magnetization just as ferromagnetics of ordinary matter does. 
This is a mechanism of ferromagnetism of the axion domain wall, which is different from a 
mechanism we have recently discussed [7] in domain walls with fermion zero modes. 

Now we show that the attractive potential for neucleons or leptons is produced by the 
axion domain wall. For the purpose we note that axion interaction || with these particles 
is specified as follow, 
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Lint = gd u a^ b ^i) + h.c. (1) 

where a is the axion field and ip is the fermion field ( for our purpose we suppose that 
the fermions are protons or electrons ). g is the coupling constant whose value depends on 
models [p,10[] of the axion. We assume that for both hadrons and leptons it is the order 



of f P Q, the breaking scale of Peccei-Quinn symmetry ( f P q = 10 10 GeV ~ 10 12 GeV ) [|J. 
When we consider a model of hadronic axion, the coupling constant with leptons is the order 
of afp^; a = 1/137. 

Suppose that the domain wall is located at x 3 = and is extending in x\ and X2 directions. 
Then the nonvanishing component of d a a is d 3 a. Dirac equation (i7 M <9 M — m + ga^^^ip = 
is reduced to 

{{E — m)a 3 — gd 3 a)u + id 3 v = 0, (— (E + m)a 3 + gd 3 a)v — id 3 u = (2) 

where E (m) is the energy ( mass ) of the fermion and we have set ip = assuming 

fermions with no momenta in X\ and x<i directions. We have used a convention in the 

Bjorken-Drell text for gamma metrices ( a 3 is Pauli metrics ). Taking account of smallness 

of gda in comparison with the fermion mass, we obtain the Schrodinger equation for the 

large component u ( v ) in the nonrelativistic limit, 

d 2 ( d 2 \ 

£U = + 9d3aa3 " >U [ £V = (~2m ~ gd?,aa3)V J ^ 

with E = m + e(E = —m — e ) where small component v = —id 3 u/2m ( u = id 3 v/2m ). 

Note that a is given by fpQ0 ( 9 is angle variable changing from to 2n around the 
domain wall ) and that the typical scale of spatial variation of 9 is given by the axion mass 
m a . Thus it turns out that the potential gd^aa^ in eq([3D for both of the fermion and the 
antifermion is attractive for the state with spin down, u = m and repulsive for the state with 
spin up, u = (Jjj . It has approximately the width of m~ l and the depth of m a . m a is about 
10 _3 eV ~ 10~ 5 eV. Hence the potential depth is quite shallow but the width is so large ( the 
order of 1 cm ) that the potential can accommodate bound states for protons and electrons. 
In such the wide potential the particles move almost freely in x 3 direction. Obviously, spins 



of these particles bounded to the wall are aligned. This causes the ferromagnetism of the 
axion domain wall. 

We may take the energies of the bound states approximately such that e = —m a + k 2 /2m 
where k% is the momentum in X3 direction. These bound states move also in x± and x<i 

— * 

directions. Then the energies of the bound states are specified with their momenta k, 

In order to calculate the magnetization we need energy spectra of these states in external 
magnetic field, B perpendicular to the wall. It is easy to derive the energies of the fermions 
and antifermions in the nonrelativistic limit, 

£M+1 ._ mo + |L +w(n+1) for£ermion 

2m 
k 2 

Ekn = — m a H — + for antifermion (5) 

2m 

where n is integer ( > ) and u is the cyclotron frequency ( uj = eB/m ). The difference in 
the energy spectra between the fermions and the antifermions comes from the difference in 
the directions of the magnetic moments of these particles bounded to the wall. 
Now we calculate a free energy of these particles bounded to the wall, 

D, = -I3- 1 lo S ( X + e^- m ~ £k '"+ l)l3 ) - /T 1 ]T log (1 + e -^ +m+£ ^ )l3 ) (6) 

where the first term represents contributions of fermions, while second one does those of 
anti-fermions. fi is the chemical potential of the baryon number or the lepton number of 
the universe. E £fc ,„ +1 <o ( E £fc ,„<o ) implies N d £ n=0 / m~ 1 dk 3 /(2n) with the condition that 
£k,n+i < ( Ek, n < ); we sum up only contributions from the bound states. ( N d is the 
degeneracy given by eBL 2 /2n where L 2 is the surface area of the wall ). 
The magnetization, M per unit volume of the wall is then given by, 

i on 

= -Turn < 7 > 

where / is the width of the wall ( / = m" 1 ). 



Let us first calculate the free energy contributed by protons and the corresponding mag- 
netization in the case of the temperature being much lower than the proton mass. Because 
we consider the case that the energy e is much lower than the temperature, we expand Q in 
terms of e~ mpl3 and e/3 and only take the leading term, 

tt=-/rv^(cosh/i/3( E + E) + sinh ^( E - E )) (8) 

where m p is proton mass. The summation is performed as follows, 

_ _ ^ r ldk 3 _ N d ly/2m p u jm a Pm£ 

= W^fJ^l + l F i/ 2 + 0(F o } _ ^ (9) 
E =^4^(^ + V/ 2 + O(F0)) U (10) 

where we have set F = m a /uj and taken the limit of the small magnetic field, F — ► oo. Note 
that a linear term in B arises only from the combination, D efe „ +1 <o — Se fe „<o> this term 
determines the magnetization. 
Thus M is given such that 

M = p- 1 smh^e- m ^ e{2mpn l a)1/2 (11) 

in the limit of the small magnetic field. Here we rewrite the chemical potential fi in terms 
of the ratio, Nb of the baryon number density to the entropy density of the universe [H; 
Nb = tib/s where % = 8e~ mpl3 smh[if3(m p /2iif3) 3j ' 2 is the baryon number density and 
s = 2n 2 g Jf /45(3 3 is the entropy density. Then it follows that 

ct^^Nb 
Mp ~ 90 (3^m p {U) 

where g* is the number of massless degrees of freedom at the temperature (3~ 1 ( = 1 ~ 100 

MeV ). 

Similarly we calculate the contribution from electrons to the magnetization of the domain 
wall. In the case the mass, m e of electron is much lower than the temperature under 
consideration. Thus we expand the free energy in terms of fi/3, m e f3 and efi, 



fi = -r 1 (io g 2-^)( e + E)-f( E -E) 

£fc,n+l<0 £fe,n<0 £fe,n+l<0 £fe,n<0 

where we have taken the leading term of the order of (3~ x and the next leading term of the 
order of \i and e. We have also taken the limit of the small magnetic field; m e is electron 
mass. Then, it follows that M e = en\/2m e m a / 4n 2 . Rewriting the chemical potential in 
terms of the ratio N L = 15/i(3/27r 2 g* of the lepton number density to the entropy density, 
we obtain 

_ ey/2m e m a g*N L 

30(3 { ' 

We find that both of the magnetization M p and M e are proportional to the baryon 
number N B and the lepton number N e , respectively. This is resulted from the fact that 
the magnetization arises owing to polarized particles and antiparticles bounded to the wall; 
their spins are aligned in the identical direction so that the sum of their magnetic moments 
are proportional to N B or N L . 

We see that the magnetization M p decreases much faster with the temperature than the 
magnetization M e does. Thus as far as we consider the temperature below lOOMeV, the 
magnetization of the axion domain wall is mainly determined by the effect of electrons if 
N B = N L . Numerically we find that 

M ~ \0'\l m " ^-Gauss (15) 

V 10- 4 eV lOOMeV lO" 10 V ; 

with g* ~ 10. This formula holds in the temperature below lOOMeV and above lMeV. ( 
When we consider the hadronic axion, magnetization due to the effect of electron bound 
states is smaller than M e by a factor of y/a. ) 

We have discussed the magnetization of the axion domain wall. The origin of the phe- 
nomena is in the axion-neucleon ( lepton ) coupling which produces the attractive potential 
around the wall. Owing to the potential the fermions are bounded to the wall with their 
spins aligned in the direction perpendicular to the wall. This dynamical mechanism has 
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been revealed by the detail calculation. But as we will show, the existence of a magnetic 
field associated with the axion domain wall can be inferred from the coupling of the axion 
and electromagnetic field, g^aE ■ B, where g 7 is a coupling constant. Let us rewrite the 
coupling in the following, 

J g~,E ■ B = J g 7 e kji djaEiA k (16) 

where we have taken the static limit, dB/dt = 0; e k ji is the antisymmetric tensor and A k is 
electromagnetic potential. Then we may identify the electric current J k associated with the 
domain wall, 

■h = 9-ytkjidjCtEi = g y e 3i kd 3 aEi (17) 

This implies a Hall effect; the electric current flows on the wall in the direction perpendicular 
to the electric field. The Hall conductivity |IT[] , a xy is given by 



v xy = 9 j J d 3 adx 3 = g J f PQ 2rm (18) 

where n is integer. By noting that g 7 = constant x fpg x e 2 /2n we find that the Hall 
conductivity is quantized; it reminds us quantum Hall effect []TTJ. This argument strongly 
suggests the existence of a magnetic field perpendicular to the wall. 

Finally as a phenomenological application we point out existence of a possible candidate 
of a primordial magnetic field fL2|,[7| generated by the domain walls. For the purpose we 
take an axion model with a heavy charged fermion and with no domain wall problem; 
the walls decay quickly. The model has axion closed strings which surround axion domain 



walls produced after QCD transition. Since the heavy fermion has zero modes |T3| on the 
strings, the strings are superconducting |TJ|]. These strings begin to carry a current Jo = Ml 
associated with magnetization M, when the walls are produced; Jo ~ 1^4 at (3 = 100 MeV 
with m a = 10 -5 e\A This current generates a magnetic field with strength Jq/Rq and with 
length of coherence Ro; Ro is the radius of the loop which is the order of the size of holizen ( 
Ro = 10 6 cm at temperature lOOMeV ). The current increases with the string loop shrinking 



in such a way that JR = JqRq because the number of the fermions ( oc JR ) on the loop 
is conserved. Then the magnetic field B = J/R becomes stronger as the radius R of the 
loop becomes smaller. We may estimate explicitely the strength of the magnetic field with 
a critical size of the coherence 10 cm jn| at the temperature 100 MeV; magnetic fields with 
smaller sizes of coherence are dissipated due to finite conductivity of the universe and do not 
remain at recombination of photons and electrons. We find that B = 10 3 Gauss with the 
size of coherence, 10 cm. This magnetic field leads to a magnetic field with strength 10~ 13 
Gauss and size of coherence 10 9 cm at the recombination ( ~ 1 eV ). This is sufficiently 
strong to be a possible candidate of a primordial magnetic field leading to galactic magnetic 
fields in the present universe. 
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